Dynamical revival of phase coherence in a many-boson system 

Shi-Jie Yang*^ and Sumin Nie^ 



''Department of Physics, Beijing Normal University, Beijing 100875, China 



o 

(N 
o 

Q 

cn 
(N 



C/2 



I 

o 



in 



X 



We study the quantum dynamics of cold Bose atoms in a double well. It is shown that self- 
trapping, as well as population oscillations are common phenomena associated to nonlinear interac- 
tions. For larger U/t, multi-particle tunneling is damped and the quantum dynamics is dominated 
by the single-particle tunneling. The many-body system can be effectively described in a truncated 
Fock space. It exhibits coherence-decoherence oscillations in the temporal evolution. We predict 
a novel phenomenon of dynamical revival and collapse of matter wave fields in optical lattices in 
regimes near the superfluid-Mott insulator phase boundary. 
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Tunneling through a barrier is a paradigm of quan- 
tum mechanics and usually takes place on a nanoscopic 
scale. The existence of long-range phase coherence, in- 
herent to the description of a Bose-Einstein condensate 
(BEC) in terms of a coherent matter wave, was demon- 
strated experimentally pHSI]- In a simple double well, 
Josephson oscillations takes place in a Bose-Einstein con- 
densate (BEC) when the initial population imbalance 
is below a critical value Above the critical value 

and as the repulsive interactions between the atoms be- 
come stronger, a phenomenon of macroscopic quantum 
self-trapping where the atoms essentially stay in one 
well is predicted 0, S] and observed for the experimen- 
tal lifetime[i,[l3|. 

As to the correlated many-body system, single-particle 
tunneling and multi-particle tunneling take effect simul- 
taneously and the collective motion of the atoms becomes 
quite complicated. Phase coherence is very difficult to es- 
tablish due to the combined contributions of single- and 
multi-particle tunneling. Virtually, self-trapping is not 
the unique phenomenon of coherent matter waves, it may 
happen in a correlated cold atom gas under similar ex- 
perimental conditions [Tl, 12 1. The physical origin of self- 
trapping is, rather than the macroscopic phase coherence 
of the system, the nonlinear interactions. Furthermore, 
population oscillations between the double wells also take 
place, given the interaction is small and the particles are 
initially imbalanced. The coherent Josephson oscillations 
in BECs is just a special case. 

In this paper, we explore the dynamical evolution of 
a correlated Bose gas in a double well by directly di- 
agonalizing the Hamiltonian in the Fock space. We at- 
tempt to clarify the roles of single- and multi-particle 
tunneling in the quantum dynamics 13 - 3] ■ It is shown 
that the so-called self-trapping, as well as the popula- 
tion oscillation, are natural consequences of the dynam- 
ical evolution. For larger U /t, the system experiences a 
coherence-decoherence oscillation. We further predict a 



phenomenon of dynamical revival and collapse of macro- 
scopic matter waves in optical lattices, where the sys- 
tem parameter is near the superfluid (SF)-Mott insulator 
(MI) phase boundary and the state is initially in the deep 
Mott regimes. This remarkable phenomenon is caused 
by the damping of many-particle tunneling whereas the 
single-particle tunneling dominating the quantum dy- 
namics. 

We concern with N bosonic atoms confined in symmet- 
rical double wells (left and right). Making the two-mode 
approximation and considering only the lowest energy 
band, the creation and annihilation operators {a\ 2 and 
ai^2 respectively) for atoms localized in either side of the 
well can be constructed. The Bose-Hubbard Hamiltonian 
is written as[l|, [l3| 
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t{a\a2 + alai) + —[hi{hi - 1) +712(^2 - 1)], (1) 



where U is the repulsion between a pair of bosons oc- 
cupying the same site, t is the hopping coefficient, and 
•hi {i — 1,2) are the number operators. Consistence of 
this model requires ^Uni{ni — 1) ^ Koi, where uji is 
the characteristic frequency of the well in the harmonic 
approximation. 

The Hamiltonian ([T]) is explicitly represented in the 
Fock basis set {|A^,0), |7V - 1, 1), • • • , |0,7V)} as 
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The general eigenstates are expressed as linear combina- 
tions of the occupation bases, \ipj) = J2k=o Cjfcl-^ — k, k) 
{j — 0,1, ■ ■ ■ , N), which correspond to the eigenvalues 
ujj. The coefficients Cjk satisfy the recursive relation 
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- V(A^-fc)(fc + l)c,(fc+i) - WiN-k + l)kcj(^k-i) + [^{N^ -2Nk-N + 2k'') - u.j\c,k = 0. (3) 
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The temporal evolution of the state is governed by the Schrodinger equation for a given initial state |5'(0)). 
I 
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The population imbalance between the two wells is 
evaluated by Z(r) = {^\{ni - ris)!*)/^ = Ef=o(l " 
2k/N)\gk{T)\'' . To depict the coherence degree of the 
system, we introduce a characteristic parameter 



{^Pjl^m and 5fe(r) 
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where Ai and A2 are the two eigenvalues of the single- 
particle density p^,.(t) = (*(r)|aj^ay |vI/(t)) {p,v = 
1,2) When a — > 1, the system is in the co- 

herent (quasicoherent) state since in this case there is 
only one large eigenvalue of matrix ppi,(r). Accord- 
ingly, a — indicates the system is in the decohcrcnt 
or fragmented state because there are two densely pop- 
ulated natural orbits. One can also use the quantum 
fluctuations of the particle number in one of the wells, 
cti(t) = y/ (\h-] (t) — (ni (T)W)/(h-] (r)), to describe the 
coherence of the svstem|21|. In the weak interaction, 
strong tunneling limit {U/t <C 1), each atom is in a coher- 
ent superposition of left-well and right-well states. The 
ground state of the system is a state with a mean number 
N/2 of atoms in each well with Poissonian fluctuations 
CTi = ^JNI2. The states in each well are quasicoherent 
states as the total number of atoms in the system is fixed. 
In the opposite limit of strong interactions or weak tun- 
neling {U jt^ 1), the tunneling term is negligible. In this 
case, the Hamiltonian is the product of number operators 
for the left and right wells. The eigenstates are products 
of Fock states {a\ = 0) and are referred as decoherent 
states. This regime is analogous to the MI phase in a 
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lattice system. 

In the following we consider two typical cases in which 
the initial state is a coherent state or a Fock state. We 
choose iV = 10 and set i = 1 as the units. 

Initial coherent states We first study the quantum dy- 
namics of an TV-particle system which is initially in a 
coherent state, 



it 
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The particles are evenly distributed in the two wells, i.e., 
Z{t) = 0, whereas the number fluctuations persist. Fig- 
ure 1 shows the temporal evolution of parameters a and 
(Ti for various interactions U/t. When U/t is small, a{T) 
fluctuates around unit and the system preserves coher- 
ence or quasicoherence. This is analogous to the SF in 
optical lattices. As U/t increases, the value of a lowers, 
indicating the coherence is gradually destroyed. How- 
ever, when U/t further increases, q:(t) again exhibits a 
nearly periodic oscillation. In the limit oi U/t ^ 1, the 
value of a oscillates between unit and zero with a period 
of T = tt/U, implying the system experiences revivals 
and collapses of the coherence. 

This phenomenon can be explained as follows. 
The initial state can be expanded as |^'(0)) ^ 
ELo'^Nia\)^~H4)''\0)- When U/t > 1, the Fock 
states (a])''|0) in each well are the eigenstates and the 
corresponding eigenenergies are Uk{k — l)/2. The sys- 
tem is then a superposition of products of the Fock states, 
which evolves independently as 
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FIG. 1: (Color online) Temporal evolution of Q(r) for the 
= 10 system initiating from the coherent state. From (a) 
to (d), U/t = 0.2,2,50,100, respectively. Periodic revivals of 
coherence is clearly enhanced for U/t^ 1. 
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FIG. 2: (Color online) Comparison of evolution of a(r) be- 
tween (a) independent and (b) full quantum dynamics at 
U/t = 100. The coherence-decoherence oscillation period is 

T = -r/U. 



There attaches a time-dependent phase factor in each 
term. When the time evolves an integer multiples of 7r/f/, 
the system recovers its initial coherent state. In between 
this period, superposition from various terms in (jS]) can- 
cels and the coherence is destroyed. Figure 2 displays 
the coincidence of the full and the independent evolving 
dynamics. 

This result has the same gradients of the experiment 
carried out by Greiner et al, where they studied the 
dynamical evolution of the matter wave field in a 3- 
dimensional optical lattice [20*1. While the system is ini- 
tially prepared in a superfluid state, the experimental pa- 
rameters are suddenly changed to the deep Mott regimes. 
A periodical collapse and revival of the matter wave fields 



were found in the temporal evolution. 

Initial Fock states We now explore our main results for 
the many-body system initially in a symmetrical Fock 
state. 



|*(0)) = |7V/2,7V/2). 



(8) 



It is found that for all values of U/t, Z(t) = 0, which 
implies the atoms are evenly distributed in the two wells. 
But the quantum fluctuations of the particle population 
in each well persist. Figure 3 shows the coherence degree 
a(T) (solid lines) and the number fluctuations (Ji{t) (dot 
lines) almost evolve synchronically. Intriguingly, with 
the enhancement of interaction, the system gradually ex- 
hibits a regular oscillating feature [Fig. 3(c), (d)]. 

This remarkable phenomenon origins from the damp- 
ing of the multi-particle tunneling due to the constraint of 
energy conservation. According to the quantum pertur- 
bation theory, the hopping coefficient for single particle 
tunneling is t. For two-particle co-tunneling it becomes 
t^ /U . As U ^ t, the single-particle tunneling process 
will dominate the quantum dynamics. With this consid- 
eration, the Hamiltonian can be effectively approximated 
in a truncated Fock space by omitting the multi-particle 
tunneling, {|iV/2-l, A^/2+1), \N/2, N/2), \N/2+l, N/2- 
1)}, where only the single-particle tunneling is relevant. 
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changed to deep SF regimes {U/t <C 1), where the multi- 
particle tunneling becomes as important as the single- 
particle tunneling, no macroscopic phase coherence will 
occur. Orzel et al 2l| had suggested the possibility of 
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FIG. 3: (Color online) Quantum dynamics for the A'^ = 10 
system with the initial Fock state |*I'(0)) = |5, 5). The solid 
lines describe a(r) and the dotted lines describe the particle 
number fluctuations o-i(r) in one of the double well. From 
(a) to (d), U/t = 0.2,2,10,80. The broken curves are a(r) 
calculated in the truncated Fock space. 



The broken curves in Fig. 3 depict the temporal evolu- 
tion of a(T) calculated with H in the truncated Fock 
space. The coincidence with the curves in the full Fock 
space (Fig. 3(c) and (d)) verify our explanation. Fur- 
ther increasing the interaction will suppress the coher- 
ence (Fig. 3(d)) as the single-particle tunneling is also 
damped. It can be proven from ^ that in the limit 
of U/t ^ 1, Q:(r) oscillates sinusoidally with a period of 
T- 2tt/U. 

Based on our result in the double well system, we pre- 
dict a novel phenomenon of dynamical revivals of macro- 
scopic matter waves for cold Bose atoms in optical lat- 
tices by comparison to the experiment by Greiner et 
al(20|. Suppose the system is initially in the deep Mott 
regimes, one suddenly changes the experimental param- 
eters to regimes near the SF-MI phase boundary, where 
only the single-particle tunneling dominates. We can ob- 
serve creation and destruction of macroscopic phase co- 
herence in the temporal evolution. At first sight, our 
prediction seems the reverse process of that of Greiners'. 
The underlying mechanism is quite different. In our case, 
macroscopic phase coherence is dynamically generated 
from a MI and the single-particle tunneling plays a key 
role. Most notably, if the system parameters are directly 



observing the temporal evolution of SF order parame- 
ter. Basically, the bosons are prepared in the number 
squeezed Mott state, and then the potential is suddenly 
reduced into the superfluid phase. The consequent evolu- 
tion of the superfluid order can be deduced from the in- 
tensity of interference patterns appearing when the atoms 
are released from the trap at sequential times. In a re- 
cent experiment, S. Will et al 22| observed multi-body in- 
teractions in the time-resolved coherent quantum phase 
revivals. 

Altman et al[2^ have studied the dynamics of bosons 
in an optical lattice using a modified coherent states path 
integral. They suggested that a system prepared in the 
unstable Mott state is expected to exhibit macroscopic 
oscillations of the superfluid order parameter. Our work 
explicitly clarified the roles of multi- and single-particle 
tunneling. We point out that the revival of coherence 
takes place only near the phase transition where the 
multi-particle tunneling is damped. If the system is di- 
rectly changed to deep superfluid regimes, then no revival 
of the coherence will be observed. This result is somehow 
contradict to one's intuition and is critical for experimen- 
tal observation. 

Finally, we briefly discuss the system evolves from an 
asymmetrical Fock state, e.g., |4'(0)) — |10,0). Figure 
4 shows the evolution of the population imbalance Z(t). 
In the weak interaction regimes (Fig.4(a)), Z(t) varies 
between —1 to 1, indicating a phenomenon analogous 
to the Josephson oscillation in BEG. We emphasis the 
population oscillations is not the Josephson effect since 
the system is not in a macroscopic coherent state. As 
the interaction increases, Z(t) remains positive. Most of 
the particles keep staying in one of the wells due to the 
constraint of energy conservation. This is the so-called 
self-trapping phenomenon. It is notable that in this case 
there is no revival of phase coherence. 

In summary, we have investigated the dynamical evo- 
lution of cold Bose atoms in a double well. We con- 
clude that self-trapping and population oscillations be- 
tween the two wells are consequences of nonlinear interac- 
tions and energy conservation instead of the macroscopic 
phase coherence. We predict a novel phenomenon of re- 
vivals and collapses of macroscopic phase coherence from 
deep MI regimes in optical lattices. Our results may be 
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FIG. 4: (Color online) Temporal evolution of the population 
imbalance Z{t) for A'' = 10 atoms starting from the Fock 
state |10, 0). From (a) to (c), the interaction U/t = 0.1, 0.5, 2. 
Self-trapping occurs when [//f > 2. 



helpful to understand how a macroscopic phase coher- 
ence can dynamically be generated from an incoherent 
state and the roles of particle tunneling in the dynamical 
evolution of quantum many-body system. 
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